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&f). Abstract 

O \ A system of equations, describing the evolution of electromagnetic 

<Z3 ■ fields, is introduced and discussed. The model is strictly related to 

^h ~*. Maxwell's equations. As a matter of fact, the Lagrangian is the same, 

O but the variations are subjected to a suitable constraint. This allows 

to enlarge the space of solutions, including for example solitary waves 
y—t ■ with compact support. In this way, without altering the physics, one 

is able to deal with vector waves as they were massless particles. The 
| main properties of the model, together with numerous exact explicit 

I ■ solutions are presented. 
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: 1 The modelling equations 

In [13], the theory of electromagnetism has been reviewed by introducing 
a suitable nonlinear set of equations, allowing for a very extended space of 
solutions, basically including any type of solitary wave. The new formula- 
tion provides a far more accurate description of wave phenomena, including 
their self-interactions. We reorganize here part of the material, providing 
additional insight. 
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We are mainly concerned with the evolution of pure electromagnetic waves 
in vacuum. Denoting by c the speed of light, in Minkowski space the model 
equations read as follows: 

— = c 2 curlB - (divE)V (1.1) 

^ = -curlE (1.2) 

E + V x B = (1.3) 

Here, E and B are the usual electric and magnetic fields, while V is a new 
velocity field. Note that there are 3 vector unknowns and 3 vector equations. 
We will require V to be of constant norm. More exactly, we impose: 

||V|| = c (1.4) 

The above condition can be eliminated in view of a more general formulation 
(see section 5). By taking the divergence of equation (1.2) we get for any t: 

divB = (1.5) 

provided the initial data are compatible with this constraint. 

By defining p = divE and by taking the divergence of equation (1.1) it is 
straightforward to get the continuity equation: 

^ = -div(pV) (1.6) 

Equation (1.6) is the first one of many other conservation laws associated 
with the above system of equations. This is quite an important prerogative, 
since, in the search of solitonic solutions, all the possible quantities must 
be preserved in the evolution. The important fact is that, even in vacuum 
p is allowed to be different from zero. Evidently, by enforcing p = 0, we 
come back with (1.1) and (1.2) to the standard Maxwell's equations, but, 
certainly, this is not the case we would like to discuss. As a matter of fact, 
it is known that solitary waves do not belong to the space of solutions of 
Maxwell's equations (we come back on this issue later in section 4). In the 
past, this led to the conclusion that quantum phenomena cannot be modelled 
by the classical equations of electromagnetism, so that, as we shall mention 
in the coming section, numerous variants of the Maxwell's system have been 
taken into account. The model we are describing here differs from the others 
for its simplicity and for many additional reasons we are going to detail. 
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2 Lagrangian analysis 



We start by introducing the potentials A and <3> such that: 



B = i 

c 



curlA 



1<9A 
c ~dt 



- V$ 



(2.1) 



In this way we guarantee equations (1.1) and (1.5). Then, we consider the 
standard Lagrangian of electromagnetism: 



It is known that Maxwell's equations are linked to the stationary points of 
the action function of C. To check this, one writes C in terms of $ and A and 
then differentiate with respect to the variations 5$ and 5 A, having compact 
support in space and time. A typical way to carry out this proof (that we 
do not report here) passes through the construction of the electromagnetic 
tensor (see for instance [17]). Thus, one discovers that (1.2) must be true 
together with the condition: p = divE = 0. 

Unfortunately, imposing p = brings to a subspace of solutions that do 
not include, for example, solitary waves with compact support. Extensions 
are then necessary and the literature is rich of results. An usual approach is 
based on a modification of the Lagrangian, by adding a further term, on the 
basis of physical considerations. With this respect, we mention the pioneering 
paper of [6] (see also [18] and [19]). For a more recent general viewpoint we 
refer for instance to [2], [3] (see also [4] for a general review). 

The procedure of adapting the Lagrangian is sometimes successful, but 
the new corresponding system of equations turns out to be heavily non- 
linear and some of the good invariance properties of Maxwell's equations 
might be lost. Theoretical results are usually addressed to the existence of 
stable solitonic solutions. Another approach is to couple some equations, 
directly derived from Maxwell's system, with other type of equations, such 
as Schrodinger or Klein-Gordon (see for example [1], [9], [11]). 

The path followed here is to recover the set of model equations (1.1)- 
(1.2)- (1.3) from the stationary points of the action function associated with 
the standard Lagrangian in (2.2). In order to avoid the condition p = 0, we 
put a constraint to the variations. In this way, since there are less degrees of 
freedom for the test functions, we come out with a larger space of solutions. 



E|| 2 -c 2 ||B|| 2 



(2.2) 
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To this end, let us impose the following relation: 

cA = $V (2.3) 

Note that the above is not a gauge condition. As a matter of fact it will 
be somehow stronger (see later). From (2.3), due to (1.4), we easily get the 
scalar constraint: 

c$ = V A (2.4) 

With the help of (2.3), we can actually obtain our set of equations. The 
formal proof of this fact is given in [13], using standard variational arguments, 
after writing C in terms of the electromagnetic tensor F^. Here we add some 
heuristic considerations. 

By standard calculus one obtains: 

C = ||E|| 2 -c 2 ||B|| 2 = E- f-I^i-W) -cB-curlA 

y c at J 

Id 1 <9E 

= — — (E • A) + - — • A - E • V$ - c A • curlB - c div(A x B) 

c dt y ' c dt y ' 

- -^ E ' A > + Ht- C H' A 

- div($E) + $ divE - c div(A x B) (2.5) 

At this point, by imposing condition (2.3) and by noting that E ■ A = 
(obtained by multiplying (1.3) by A), the above equation becomes: 

C = ||E|| 2 -c 2 ||B|| 2 = -div($(E + VxB 

1 / (9E \ 
+ - I — - c 2 curlB + (divE)Vj - A 



(2.6) 



We now observe that £ vanishes when: 

||E|| = c||B|| (2.7) 

Consequently, if (1.3) and (1.1) hold true, we are exactly in the situation 
C = 0, that is: all our waves will have the norm of E equal to that of cB. 
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More technical (but not too difficult) is to show that C = is actually a 
stationary point of the action function associated with C 

Finally, we remark that the constraint (2.4) might be included in the 
Lagrangian, as a penalty term, in the following way: 

C = IIEII 2 -c 2 ||B|| 2 + -(A-V-c$) (2.8) 

c 

The added term is similar to the one we would have in presence of moving 
charged particles in an electromagnetic field (see, e.g.: [17]). 



3 Eikonal equation 

With a simple analysis, more general properties of the solutions can be re- 
covered. For example, by expressing E and B in term of the potentials, one 
can prove that: 



d 

c 2 (E + V x B) = ~^(cA - $V) + V div(cA - $V) 

+ V x curl(cA - $V) - cV • ^ + ^ V|| V|| 2 (3.1) 

where DA/Dt = dA/dt + (V ■ V) A is the substantial derivative of A along 
the velocity vector field V. Thanks to (2.3) and (1.4), one concludes that: 

DA 

E + VxB = <^ V-— - = (3.2) 
Moreover, the following relations hold: 

E-V = EB = VB = V = c „ E X B „ (3.3) 

||ExB|| v ; 

The first one is obtained by scalarly multiplying (1.3) by V. The second one 
is similarly obtained by multiplying (1.3) by B. The third and the fourth 
ones are consequence of the two previous orthogonality relations together 
with conditions (1.4) and (2.7). Thus (E, B, V) is a right-handed orthogonal 
triplet. 
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Thanks to the above relations one can prove a classical results regarding 
the time derivative of the energy: 

1 

- — (||E|| 2 + c 2 ||B|| 2 ) = c 2 (curlB ■ E — curlE • B) = — c 2 div(E x B) (3.4) 

where E x B is the Poynting vector. 

Additional conclusions are obtained by imposing conditions on the field 
V. For example, one may require that V is not subjected to transversal 
acceleration: 

DV dV 

_ = _ + (V .V)V = (3.5) 

A way to satisfy (3.5) is to suppose that V is an irrotational stationary 
field. This means that V = for some scalar potential In this way: 
DV/Dt = V(d^/dt + ±V||^|| 2 ) = 0, yielding (3.5). With this choice, the 
most important consequence is the following eikonal equation: 

||W|| = c (3.6) 

directly obtainable from (1.4). This ensures that our solutions develop ac- 
cording to the rules of geometrical optics. Note also that, if we use (2.3) 
and (3.5), one gets: cDA/Dt = D{$V)/Dt = V(D<Z>/Dt). Therefore, the 
condition at the right-hand side of (3.2) can be replaced by: 

_ = _ + V .V* = (37, 

In conclusion, when V is irrotational and DQ/Dt = 0, we automatically 
have (1.3) and (3.6), showing that the wave-fronts evolve as prescribed by 
geometrical optics. 

If, in addition to these hypotheses, one also requires the following Lorenz 
gauge conditions on the potentials: 

1 <9$ 

divA + -— = (3.8) 
c at 

from (2.3), we also get another continuity equation: 

— = -div($V) (3.9) 



which is stronger than demanding DQ/Dt = 0. 

Our final purpose is to show that it is possible to explicitly compute infi- 
nite interesting solutions of the system (1.1)-(1.2)-(1.3), realizing an extended 
range of conservation laws. 

4 Explicit solutions 

In the cartesian reference frame (x,y,z), let us take the following solution 
candidates: 

E= (E 1 (x,y)g(z-ct), E 2 (x,y)g(z - ct), o) 

B= (B 1 (x,y)g(z-ct), B 2 (x,y)g(z - ct), o) (4.1) 

representing two field distributions, modulated by the function g, laying on 
parallel planes and shifting in the direction of the z- axis at the speed of light. 
The functions E\, E 2 , B\ and B 2 are smooth on the whole plane (x,y) (they 
may allowed, for example, to be zero outside a domain of finite measure). 
The function g is also smooth. 

By setting V = (0, 0, c) and by direct substitution into the equations 
(1.1)-(1.2)-(1.3), the vector fields shown in (4.1) are solutions under the fol- 
lowing assumptions: 

E x = cB 2 E 2 = -cB 1 divB = dB 1 /dx + dB 2 /dy = (4.2) 

for any choice of g. With this setting, one can actually find a function A 3 
and construct the two potentials: 

$ = A 3 (x, y)g{z - ct) A = (o, 0, A 3 (x, y)g{z - ct)) (4.3) 

It is straightforward to check that all the conditions (2.3), (2.7), (3.3), (3.4), 
(3.5), (3.6), (3.7), (3.8), (3.9) are verified. Thus, depending on the arbitrary 
functions A 3 and g, we can build infinite solutions of our system. If A 3 
and g have compact support, we get solitary electromagnetic waves shifting 
unperturbed at the speed of light. Note that p = divE ^ and (1.6) becomes 
the trivial transport equation: 



If we instead start from the general setting in (4.1), and try to solve 
the full set of Maxwell's equations (i.e., including the additional condition 
divE = 0), there are no chances of getting interesting nontrivial solutions. In 
this case we can forget about the velocity vector V. By direct substitution, 
one easily deduces that both the functions E\ — %E 2 and B 1 — %B 2 , where 
% is the imaginary unit, must be holomorphic (entire) on the whole complex 
plane x + iy (see, for instance [10]). By the Liouville's theorem: if an entire 
holomorphic function is bounded, then it is a constant, we entail that there 
are no bounded continuous electromagnetic fields of the form (4.1), having 
finite energy and solving the whole set of Maxwell's equations, with the 
exception of E = and B = (note that classical plane waves have not 
finite energy). We recall that solutions may exist if we assume that E and B 
do not belong to the tangent plane of the advancing front. In this situation, 
however, the Poynting vector E x B, indicating the direction of the energy 
flow (see (3.4)), is not lined up with the direction of movement, in contrast 
with the rules of geometrical optics. This brings to a diffusive behavior. 

We can now transform the equations (1.1)-(1.2)-(1.3) in spherical coordi- 
nates (r, 6, 0) and take the following fields, distributed on the tangent planes 
of spherical wave-fronts: 

E= 1(0, E 2 {6,<t>)g{r-ct), E 3 (9,cf>)g(r - ct)) 

B = ^(0, B 2 {6,<t>)g{r - ct), B 3 (9,<f>)g(r - ct)) (4.5) 

With this choice, the energy density ||E|| 2 + c 2 ||B|| 2 remains constant when 
integrated over any spherical surface. Similarly to the previous case, by 
defining V = (c, 0, 0), we can get infinite solutions, provided they satisfy: 

E 3 = -cB 2 E 2 = cB 3 divB = d(B 2 sin B)/d0 + dB 3 /d(f) (4.6) 

allowing us to build the potentials: 

$ = A 1 (9,<f>)g(r -ct) A = (A 1 (9,cf>)g(r - ct), 0, o) (4.7) 

where A\ and g are arbitrary. As before all the conditions (2.3), (2.7), (3.3), 
(3.4), (3.5), (3.6), (3.7), (3.8), (3.9) are verified. By taking A^B, <p) = -cos9 
and g(s) = sinus, among the solutions we recognize the one, corresponding 
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to the monochromatic field generated by an infinitesimal dipole, usually em- 
ployed in applications (see [5], [7], [15]): 

„ / sin^ . \ „ / sin# . ,\ 

E= [0, —-smcu(r - ct), 0) B = (0, 0, — ^— sin u(r - ct)j (4.8) 

Searching instead for solutions of the Maxwell's system, one finds out 
that the complex functions Ei — 1E3 sin 9 and B2 — 1B3 sin 9 should be 
holomorphic on the Riemann sphere. Hence, E and B are bounded if and 
only if they are zero. Therefore, the only possible solutions of Maxwell's 
equations in vacuum, having bounded electromagnetic fields laying on tan- 
gent planes of spherical fronts are identically zero. As a consequence the one 
in (4.8) cannot be solution (in fact divE 7^ 0). The same arguments can 
be applied to any closed, bounded, compact, oriented surface. Non trivial 
solutions, for example the celebrated Hertz solutions (see [16], [7], [15]), are 
instead possible by allowing the fields to have radial components E\ and B ± . 
In this way, however, the Poynting vector E x B is not aligned with the 
direction of movement and the evolution of the corresponding wave-fronts do 
not comply with the rules of geometrical optics (i.e.: the surfaces obtained 
as the envelope of the electromagnetic vectors, representing the wave-fronts, 
are not spherical), so that many of the conservation laws reported here are 
not satisfied. 



5 Concluding remarks 

With the new model equations we have an extended range of solutions, not 
available in the Maxwellian case. In particular, this includes wave-packets of 
almost any form, both from the viewpoint of the shape of the wave-fronts and 
the information written on them. These wave-packets travel unperturbed at 
speed c, along the direction of the vector V, having constant norm. If V 
is stationary and irrotational, then the eikonal equation is verified and the 
wave-fronts perfectly follow the laws of geometrical optics. The equations 
are compatible with the stationary points of the standard Lagrangian of the 
electromagnetism. In this new framework, the differentiation has to be taken 
with respect to potentials subject to a certain constraint, that, thanks (3.1), 
is equivalent in practice to impose (1.3). 
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At this point, something has to be said about the physical implications. 
First of all, we recall that in [13] equation (1.3) is generalized as follows: 



DY 



/i(E + V x B) 



Wp 



(5.1) 



Dt 



P 



which is the Euler equation for compressible fluids (recall the continuity 
equation (1.6)), with an electromagnetic type forcing term. The scalar p 
is a suitable pressure and p is a constant whose dimension is charge/mass. 
This approach combines the evolution of electromagnetic entities with that 
of a (non material) inviscid fluid. Some numerical experiments, regarding 
the interaction of waves with matter (i.e., the diffraction due to the passage 
of a photon through a small hole) are examined in [14]. Other numerical 
experiments, concerning electromagnetic waves trapped in bounded regions 
of space (basically, forming vortex rings) are reported in [8]. In this new 
context, that is better suited for a general relativity framework, we must 
drop condition (1.4) and replace it with an eikonal equation in a suitable 
metric space. In order to maintain the exposition at basic level, we do not 
add further considerations concerning this generalization. 

So far, we studied the special case when p = and DY / Dt = 0. In this 
circumstance, the fluid moves unperturbed. From the electromagnetic point 
of view, equation (1.3) says that the "Lorentz force" acting on the wave is 
null. This is in agreement with the fact that our solitons freely travel, without 
the influence of external factors. 

Equation (1.1) is the Ampere law, without explicit external currents, 
but with a sort of electric density developing together with the wave and 
satisfying automatically the continuity equation (1.6). The presence of a 
nonvanishing p is not in contrast with the Gauss's divergence theorem. If 
our travelling soliton has a bounded support included in a solid region Q, 
then the integral J sn E • n = J n p is equal to zero. Therefore, seen as a 
whole, our soliton is not a charge (it does not even emit electromagnetic 
fields during its movement), although, inside, there are points where p ^ 0. 

In order to deal with equation (5.1) we may generalize the Lagrangian in 
(2.8) by setting: 



C = ||E|| 2 -c 2 ||B|| 2 + £(A-V-c$) + — Vc 2 - ||V|| 2 - V - (5.2) 
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We come back to the special case of equations (2.8) by imposing p = and 
recalling (1.4). This approach takes also into account relativistic effects. We 
did not carry out any analysis concerning the above Lagrangian. Perhaps, 
the techniques introduced in [4] to study the general Lagrangian £ U! s,a,<s> m &y 
also be applied in our situation. 

To conclude we say that the most important achievement, realized with 
the new model equations, is to be able to produce, maintaining a physical 
meaning, vector waves that display all the peculiarities of standard electro- 
magnetic emissions. In addition, waves with the characteristics of a photon 
are allowed by the model, providing the link between the evolution of elec- 
tromagnetic entities and that of classical mechanical bodies. Finally, let us 
point out that the new system of equations can be proved to be invariant 
under Lorentz transformations (see [13], section 2.6), can be written in co- 
variant form (see [13], chapter 4) and is naturally linked to the derivative of 
the electromagnetic stress tensor (see [13], section 4.2, and [12]). 
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